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Abstract. We present an approach to quantum dynamical lower bounds for 
discrete one-dimensional Schrodinger operators which is based on power-law 
bounds on transfer matrices. It suffices to have such bounds for a nonempty 
set of energies. We apply this result to various models, including the Fibonacci 
Hamiltonian. 



Consider a self-adjoint operator H on a separable Hilbert space Ti — £ 2 (Z d ) or 
£ 2 (N). The dynamical evolution of an initial state ip E Ti, = 1, is given by 
tp{t) = exp(—itH)ip. We shall denote ip(t,n) — (ip(t),5 n ), where B — {S n } is the 
canonical basis of Ti labelled by n £ 7L d or n € N. One usually takes ip so that 
■0(0, n) = ip(n) be well localized (fast decaying at infinity). For example, one can 
take ip = S± in the one-dimensional case. While being localized at t = 0, the wave 
packet in general spreads with time over the basis B. It is convenient to consider 
the time-averaged quantities 



(1) a(n,T) = - \ip(t,n)\ 2 dt, ^ a(n, T) = \\4'\\ 2 = 1 for all T > 0. 



There exist basically two possibilities to characterize the spreading of the wave 
packet. First, one can consider the upper and lower rates associated to the fastest 
(or the slowest) part of the wave packet. Let 



France 



1. Introduction 




n 



a(n,T), 7 >0. 



n:\n\>Ti-2 



For the fastest part, one defines 
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(2) 7 -= S up^ 7 >0| Hmmf^|^ = 

(3) 7 + = sup i 7 > | limsup = 



(with convention logO = — oo). Assume that 7 ± > 0. Then it follows from this 
definition that for any < v < r\ > 0, 



(4) a(n,T)>CT- r > 

n:|n|>TT"-" 

for all T > 1 with a uniform (in T) constant C. For 7 + , a similar bound holds on 
some sequence of times Tk — > +oo. One could take a slightly different definition of 
7 ± so that one has for any < v < 7", 

(5) £ a(n,T) > > 

r^lnl^TT - -" 

(and similarly for 7 + ), but the meaning of numbers 7^ will be essentially the same. 

The spreading rates for the slowest part of the wave packet can be defined with 
summation over {n : \n\ < T 7 } and taking inf{7} in (||), (|||) (see also Q for a 
slightly different definition). In the present paper we shall be interested, however, 
only in the fastest part of the wave packet. 

To determine the numbers 7 ± for concrete quantum systems, one should be able 
to relate it to the spectral properties of H . What is in fact known are some general 
lower bounds for 7 ± . Although not stated in this form, it follows from the proofs 
of §, H that 



(G) 

and from the proof of p0| that 



7 > -dinifj^), 



(7) 7+ > -dimpOty) 

(with d > 1 in the case of l 2 (Z d ) and d = 1 in the case of Z 2 (N)). Here ^ is the spec- 
tral measure associated to the state ip and the operator H, and dim# (//), dimp(/i) 
denote the (upper) Hausdorff and packing dimensions of the measure /i, respectively. 
For different definitions and relations between them, see, for example, j|, [Tt], [5?J. 
In particular, the Hausdorff dimension is determined by the most continuous part 
of the measure 

dim#(/i) = /i-ess-sup \~(E), 

where A - (E) is the lower local exponent of /i, 

WF x ,. • f logy{[E-e,E + e}) 

A (E) = hmmf , E G suppu. 

e^o loge 

For the packing dimension, we have 

dimp(/i) = /i-ess-sup X + (E), 
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where X + (E) is the upper local exponent of the measure, 

x + log/x([-B — e,E + e]) _ 

A + (£) = hmsup : 3 -, E e suppu. 

e^O l0g£ 

One can observe that for Hausdorff dimension, the result is slightly stronger than 
(||), namely, (ph holds with 7~ = dimij(/i^,). 

If one knows the continuity properties of the spectral measure (encoded in 
X~(E)) and one has also some information about the decay of the generalized 
eigenfunctions (n, E), associated to the state tp, then one can improve the lower 
bound (§) for 7- ||, ||. 

More involved quantities which describe the fastest portion of the wave packet 
are the time-averaged moments of the position operator: 



(8) (\X%)(T) = ± f Ty £\n\*m,n)\ a = $>|*o(»,r)> p>0 

1 J n 

(where |n| is the norm in 7L d in the case of l 2 (1 d )). One can define the associated 
upper and lower growth exponents, 



(3 (p) = liminf 



log(|XQ(T) 



t^+oo logT 

and in a similar manner for f3 + (p). The numbers f3 ± (p) depend, in general, on the 
state ip, but we will leave the dependence implicit. It is clear that (jj) implies 

(\X\^)(T) > CT"^"-")-" 

for any v > 0, 77 > uniformly in T > 1 and thus 0~(p) > P7~ ■ Together with (|^) 
this yields 



(9) /T(p) > ^dim H {^), 

the bound most often used to bound from below the moments of the position 
operator. Similarly, one always has (3 + {p) > p^ + > §dirnp(/^,). 

It is important to observe that strict inequalities 0~(p) > pq~~ , (3 + (p) > p^ + 
may occur. This is possible if the wave packet has polynomially decaying tails at 
infinity. Assume, for example, that for some 7 > 7" and some r > 0, 



(10) a(n,T)>CT~ T 

n:\n\>T-i 

uniformly in T > 1. Then f3~ (p) > pj — r. This bound is better than {3~(p) > P7 _ 
for p large enough. 

General lower bounds for the time-averaged moments which take into account 
(in a somewhat hidden form) polynomial tails, are obtained in [|| and The 
proofs are based on the spectral theorem and develop the ideas of Guarneri 
|l9| . The obtained lower bounds are expressed in terms of spectral measure 
H> 13 01 m terms of both spectral measure and generalized eigenfunctions 
u^(n,E) jj7|. To apply them to concrete quantum systems, one should have a 
good knowledge of n^([E — e, E + e]), EG supp/j^, (and also of the functions 
Sn(E) = J2 n -\n\<N \ u 4>( n iE)\ 2 in the case of the mixed lower bounds from Q). 
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Such a kind of information is difficult to obtain in the multidimensional case. How- 
ever, in one dimension for operators H of the form 

(11) Hi/}(n) = ip(n - 1) + ip(n + 1) + V(n)ip(n), 

there exist rather effective methods based on the study of solutions to the formal 
eigenfunction equation Hu = Eu, E G K. In particular, the growth properties of 
the transfer matrix T{n, E) associated with this equation are closely related with 
the spectral properties of operator H. For simplicity, from now on we shall always 
take ip = 8\ and we write simply fi instead of ns 1 ■ 

Our purpose is not to give a detailed list of results obtained in this field. We will 
mainly be interested in models with power-law growth of the norm of the transfer 
matrix. Roughly speaking, the slower the growth of ||T(n, E)\\, for E E A and 
as \n\ — > oo, the more continuous is the restriction of u on A. What will be of 
particular interest for us in the present paper, are the two following results. First, 
if \\T(n,E)\\ < C(E)\n\ 7 > for /i-a.e. E G A, fi(A) > with r/ G [0,1/2) and finite 
C(E), then @ || 

H{[E -e,E + e})< D{E)e 1 - 2, \ E e A, 
so that the measure is (1 — 2?7)-continuous on A. In particular, 

(12) diiM/i) > 1 - 2 V , pr(p)>p(l-2r}). 

More generally, assume that every solution of the equation Hu = Eu, E G A 
obeys 

(13) diE)^ 1 < Hn,E)\ 2 <C 2 (E)L 2q2 

n:\n\<L 

with positive finite C\(E),C2(E) and < q± < qi < +oo. Then the measure is 
291/(91 + 92) -continuous on A In particular, dim#(^) > 2qi/(qi + q 2 ). 

One can observe that the condition \\T(n, E)\\ < C(E)\n\ a with some a > implies 
92 < a+ 1/2. 

The polynomial upper bound on the norm of the transfer matrix also allows 
one to bound from below fJ,([E — e,E + e]). Such a bound can be used, following 
general results of ^| , to obtain non-trivial dynamical bounds for the moments 
(/3~(p) > 0) for some systems with pure point spectrum [jl8). Since in this case 
dim^f(/i) = 0, these bounds cannot be obtained by usual methods based on the 
continuity properties of the measure. (Another example of this kind, studied by 
Jitomirskaya et al. in ^5|, will be discussed later on). 

As was said above, the general results of || ^7) were obtained using the spectral 
theorem for the operator H and thus by representing all the quantities of interest as 
some integrals over the spectral measure. In the present paper we propose another 
method based on integrals over Lebesgue measure rather than over /i^,. One can 
bound from below rather directly the sums Yln-\n\>L a ( n > (and thus the time- 
averaged moments (\X\ P )(T)) using the Parseval equality (32). Namely, for any 

1 r + °° 1 f 

(14) - ^ e - 2 </ T |^,n)| 2 di=— Jj(R(E + i /T)i J )(n)\ 2 dE. 
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Here, n £ 7L d or n S N, and R(z) = (H — zl)^ 1 is the resolvent of H. The 
integral on the l.h.s. in ( [l4| ) is very close to the Cesaro time-averaged quantity 
a(n,T). Therefore we modify the definition of a(n,T) in ([[]) and of (\X\ P )(T) in 
(||) accordingly (this is done in this way by many authors) , 

(15) a(n,T) = - / e^^n)] 2 dt, (\X\^)(T) = ^ \n\ p a{n,T), 

and work with these definitions in what follows. One can easily see that this does 
not change the numbers (3^(p), provided that the moments |X|^(i) do not grow 
faster than polynomially (which is true in most applications and, in particular, in 
all applications considered in the present paper). 

To obtain a non-trivial dynamical lower bound, it is sufficient to bound from 
below Xm-|n|>n \R{E + i/T)t[;(n)\ 2 for some 7 > for E from some set of positive 
Lebesgue measure. In our one-dimensional case with ip = 8%, the resolvent R(E + 
i/T)5i{n) is closely related to the corresponding transfer matrix T(n, l;E + i/T), 
which is in turn close to T(n, l;E) if T is large and n is not too large. Roughly 
speaking, if ||T(n, 1; E)\\ is bounded from above polynomially in n for some E, then 
R{E + i/T)8\{n) decays no faster than polynomially for n not too large (depending 
on T). The following result will be proved in Section 2: 

Theorem 1. Let the operator H in £ 2 (Z) or l 2 (N) be given by (|Tl| ) (with Dirichlet 
boundary conditions at in the half-line case). Suppose that for some K > 0, 
C > 0, a > 0, the following condition holds: 

For any N > large enough, there exists a nonempty Borel set A(N) C K such 
that A(N) C [-K,K] and 

(16) \\T(n,m;E)\\ < CN a VE € A(N), V n,m : \n\ < N,\m\ < N 

(resp., with 1 < n < N, 1 < m < N in the case of £ 2 (N)). 

Let N(T) = T 1 / (1+Q ^ and let B(T) be the 1 /T '-neighborhood of the set A(N(T)): 

B(T) = {EeR:3E' £ A(N(T)), \E - E'\ < l/T). 

Then for all T > 1 large enough, the following bound holds: 

(17) Yl a(n,T)>~\B(T)\N l -* a (T), 

\n\>N(T)/2 

where C is some uniform positive constant and \B\ denotes the Lebesgue measure 
of the set B. 

In particular, for any p > 0, one has the following bound for the time-averaged 
moments of the position operator: 

(18) (\X\l)(T)>^\B(T)\N^- 2 «(T), 
where (|X|^)(T) is defined as in (fL5|). 

In the statement of Theorem |l| it is important that the constant C in (16) 
be independent of TV and of E € A(N). If ^4(^V) consists of a single point Eq 
(independently of N) , the statement of the theorem is still non-trivial. Namely, 
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Corollary 1.1 (One-Energy Theorem). If 

(19) \\T(n,m;E )\\<C(E )(\n\ + \m\r 

for some Eq, uniformly in n,m, then 

„ , . p — 1 — 4a 

0-(p)> ~ . . • 
1 + a 

If one only knows that \\T(n, 1; Bo)|| < C(Eo)\n\ n , then one can take a = 2r/, so 
that 

(20, 

It is interesting to compare this bound with (|l2"|), stating for < r\ < 1/2, that 
P~{p) > p(l — 2??). One observes that J2C| ) yields a better result for p large enough 
(and one needs the polynomial upper bound for the transfer matrix only for a single 
point and not on some set of positive measure). Moreover, in the case rj > 1/2, 
the bound (|2(j) is always non-trivial for p large enough, while (f|) gives nothing. 
In this case the spectrum of H may be pure point with polynomially decaying 
eigenfunctions. 

Another interesting consequence of Theorem [l] is the following: 
Corollary 1.2. If there exist C > 0, fi £l, 6 > 1 such that 

\\T(n,l;E)\\<C 
for all n,E such that \n\ ■ \E — Eq\ 6 < 1, then 

rip) >p- i/o. 



Taking 9 = 2 in Corollary 1.2, we obtain the main dynamical result of ]25| |. 

Finally, we note that all the results mentioned above are stable with respect to 
finitely supported perturbations. This is of interest since all the previous dynamical 
lower bounds were based on dimensionality properties of spectral measures which 
are very sensitive to such perturbations. 

Corollary 1.3. Suppose the potential Vq is such that the operator Hq with potential 
Vq satisfies the hypothesis of Theorem Let W : Z — > M (resp., W : N — ► M) be 
a finitely supported perturbation and H = Hq + W . Then the operator H satisfies 
(0) and ©. The same kind of stability holds for Corollaries \l.f\ and 



We shall apply the results above to various models from one-dimensional qua- 
sicrystal theory (cf. Jl^]): the Fibonacci Hamiltonian, the period doubling Hamil- 
tonian, and the Thue-Morse Hamiltonian. More applications will be discussed in a 
forthcoming publication . 

The Fibonacci potential V : Z — > R is given by 

(21) V(n) = \x[i- u ,i)(nu modi). 

Here, A > and ui is the inverse of the golden mean, that is, ui = (V5 — l)/2. For 
a given A, let 



(22) 



C x = 2 + V8 + A 2 . 
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We shall prove 
Theorem 2. (a) For every X,p, we have 



0~(P) > 



p — 1 — Act 



a 



with 



(23) 



a(A) 



21og(C A (2C A + l) 2 ) 



logo; -1 

(b) For every p and every A > 4, we have 

p — 7 — 3a 



0-(p) > 



1 



a 



wh 



ere a is as in 



|) and 



7 = 7(A) 



log(2A + 22) 



1 < 1 



(24) , 

logcj 1 

Let us compare the result of Theorem with previously known lower bounds 
for the moments. In fact, all these bounds are based on (H). The lower bounds for 
the Hausdorff dimension for the Fibonacci Hamiltonian [24 , 26 in turn have all 
been proved using (p~3[) . For example, in p6f it was proved that 



(25) dim ff (X) > 2n/(n + a + 1/2), 

where 

log(VTr/4) _ m9fi 

« = Tl — ? -n ~ 0.0126 
51og(w !) 

and a is given by (23[). For large p, the dominant expression in the lower bound for 
(3~ (p) in Theorem 3is p/(l + a) which is clearly better than 2np/(n + a + 1/2). We 
therefore see that our result improves the previously known lower bound for large 
P- 

We also note that Corollary 1.3 is particularly interesting in this context. The 
Fibonacci Hamiltonian has spectrum a of zero Lebesgue measure |fl [36| and hence, 
by a suitable application of the Simon- Wolff argument p4| , a generic rank two 
perturbation (at two consecutive sites) will produce a spectral measure which is 
entirely supported away from a. More precisely, the spectral measure of the per- 
turbed operator will be supported on a countable set of eigenvalues with a being the 
set of accumulation points of these eigenvalues. In particular, a generic rank two 
perturbation may turn the spectral measure of the Fibonacci Hamiltonian, which 
is known to be a-continuous for some a > 0, into a pure point measure. In this 
situation, previous bounds do not give anything for the perturbed operator, while 
our bound for f3~ (p) is stable with respect to such a perturbation. 

Our next application is to the period doubling Hamiltonian which was consid- 
ered, for example, in |6[ ||, O [2lJ]. Let fl p d be the (two-sided) subshift generated 
by the substitution i— > 01, 1 i— > 00 and define, for some u> £ f2 p d and a given cou- 
pling constant A, the potential V\ yUJ by V\^(n) — \ui n . Is is known that for every 
A / and every u> S fi p d, the operator with potential Va,w has purely singular 
continuous spectrum (see [[l3| for this result and || ||, O ^l| for earlier partial 
results). We shall show the following: 
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Theorem 3. For every A and every lu £ Q p d, we have 

r(p)>^. 

This result is of interest for various reasons. First of all, this is the first dynamical 
result for this model that goes beyond the RAGE theorem (i.e., the dynamical 
result that follows from purely continuous spectrum). Indeed, apart from singular 
continuity, nothing on either dimensionality or dynamics was known rigorously. The 
difference to the Fibonacci case, where, as discussed above, earlier results in these 
directions were known, stems from the absence of an invariant for the trace map. 
More concretely, it could not be shown that the trace map orbits remain boun ded 
for (sufficiently many) energies from the spectrum. With a result like Corollary [O] 
at our disposal, very weak results on the period doubling trace map suffice already 
to allow us to establish the result above. Secondly, observe that the dynamical 
bound we obtain in Theorem || is independent of the coupling constant. This is in 
sharp contrast to the Fibonacci case where our bound (and previous ones) becomes 
worse as the coupling constant is increased. 

Finally, we consider the Thue-Morsc Hamiltonian which was studied, for ex- 
ample, in [Q, ||, |l6[ |2lj. Let fl tm be the (two-sided) subshift generated by the 
substitution i— ► 01, 1 i— * 10 and define as above, for some u € f^tm and a given 
coupling constant A, the potential V\ tU} by V\ M {n) = Xui n . We shall show the 
following: 

Theorem 4. For every A and every lo G fitm? we have 

p-(p)>p-l. 

As in the previous theorem, the bound is A-independent and the first dynamical 
result for this model. Furthermore, in this case the spectral type has not even been 
identified in all cases. It is known that for every A ^ and every u> € ^tm, the 
absolutely continuous spectrum is empty (this follows from results of Kotani [p8| 
and Last and Simon |30|l ) , but absence of eigenvalues is only known generically that 
is, for a dense Gs set of <x> € S] tm (this was shown by Delyon and Peyriere |lq] and 
Hof et al. (2lJ], using different methods). This illustrates the theme of this paper 
in a nice way: While we cannot say anything about the dimensionality of spectral 
measures, we can nevertheless prove very strong dynamical bounds. 

The organization of the article is as follows: Theorem [l] and Corollaries LI L3 



will be proved in Section 2. In Sections 3-5 we shall then apply these criteria 
to the Fibonacci model, the period doubling model, and the Thue-Morse model, 
respectively. 

2. A New Criterion for Quantum Dynamical Lower Bounds 

Let H — —A + V be a self-adjoint operator on £ 2 (Z) or on £ 2 (N) (with Dirichlct 
boundary condition in the latter case). Given some real function V(n), it is formally 
defined by 

Hip(n) = ip(n - 1) + ip(n + 1) + V(n)ip(n), n e Z 

and by 

Hip(n) = ip(n - 1) + ip(n + 1) + V(n)ip(n), n > 1, 
#V(1)=V>(2)+IW(1) 
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respectively. Let z be some complex or real number. We define the transfer matrix 
associated to the operator H as follows: 

{A(n, z)A(n — 1, z) ■ ... ■ A(m + 1, z) if n > to, 

I if n = rn, 

T _1 (m, n; z) if n < to, 

where 

A(n, 2 )=( 2 -f("> "J 

For 2 £ C with Im z ^ 0, define 

<j> = R(z)6i = (H - zl)- 1 5i. 
By the definition of the resolvent of the operator H, we have 

(27) 4>(n-l) + (f>(n + l) + (y(n)-z)<f>(n) = 0, n ^ 1, 

(28) 0(O) + #2) + (^(l)-*Ml) = 1 
in the case of £ 2 (Z) and 

(29) <£(n - 1) + <f>{n + 1) + (V(n) - z)0(n) = 0, n > 1, 

(30) 0(2) + (V(l) - z)<l>(l) = 1 

in the case of £ 2 (N). Consider the vectors Q(n) — (<fi(n+ 1), 4>(n)) T . One can easily 
see that $(n) = A(n, z)$[n — 1) for n ^ 1. Therefore, 

(31) $(n) = T(n,l;z)$(l), n > 1 

in both cases. Moreover, in the case of £ 2 {1>), one has the identity 

(32) $(n) = T(n, 0; z)$(0), n < 0. 

Since det A(n, z) = 1, the same is true for T(n,m; z), so that ||T -1 (n, m; z)\\ = 
\\T(n,m; z)\\ for all n,m,z. Therefore, (|3l|) and (B2|) imply 



(33) ||o>(n)||>||T(n,i;z)ini<i>(i)||, n>l 

in both cases and 



(34) ||$(n)||>||T(n,0;«)||- 1 ||$(0)||, n < 0, 

in the case of £ 2 (Z). It is clear that to get a lower bound for the resolvent, we 
need an upper bound for the norm of the transfer matrix for complex values of z. 
Usually, one has such bounds only for real values of z, so we should establish 
relations between these two cases. 

Lemma 2.1. Let E eR,N > 0. Define 

K(N)= sup \\T(n,m;E)\\ 

\n\<N,\m\<N 
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in the case of I 2 [J*) and 



L(N)= sup \\T(n,m;E)\\ 

Kn.m<N 



B = 



in the case of ^ 2 (N). Let 6 G C. The following bounds hold: 

(35) \\T(n,l;E+5)\\<K(N)exp(K(N)\n\\5\), l<n<N, 

(36) \\T(n,0;E + S)\\<K(N)exp(K{N)\n\\6\), -N < n < 0, 
in the case of 1 2 {X) and 

(37) \\T{n, l;E + 8)\\ < L(N) exp(L{N)\n\\6\), l<n<N, 
in the case of £ 2 (N). 

Proof. The proof is virtually identical to the proof of Theorem 2 J in |33|]. For 
example, for any n with 2 < n < N, one can write the identity 

n-l 

T(n, l;E + 5) = T(n, 1;E) + S^2 T{n, j + 1; E)BT(j, 1;E+ 6), 

i=i 

where 

' 1 


By iteration, using the fact that \\T(n,m; E)\\ < K(N) for all 1 < n, m < N, one 
can show (|35|). The same proof yields (|3?]). The bound ( ^6|) can be proved in a 
similar manner. □ 

Proof of Theorem [|. We give the proof in the case of ^ 2 (Z), the proof in the case 
of £ 2 (N) is analogous, with some aspects being even simpler than in the whole-line 
case. 

The basic idea is the same as in [p8| in the case of a model with a sparse potential. 
First, using Parseval's identity [B2L one can write for any n, 

If 00 e f 

(38) -/ \exp>{-itH)5 1 {n)\ 2 exp(-2t/T)dt= — / \R(E + ie)S 1 (n)\ 2 dE, 
T J 2n J & 

where e = 1/T. For given E, e, we shall write 4>(n) = R(E + ie)S\(ri). The bounds 
({H and (||) allow us to bound ||$(n)|| 2 = \(p(n + 1)| 2 + \(j)(n)\ 2 from below if we 
find upper bounds for \\T(n, 1; E + ie)\\ or \\T(n, 0; E + ie)\\ , provided that ||$(1)|| 
or ||$(0)||, respectively, is not too small. 

Let T > 1. Define N = N(T) = T 1 /( 1+Q ). The assumption of the theorem gives 

K(N)= sup \\T(n,m;E')\\ < CN a , 

\n\<N,\m\<N 

provided that E' £ A(N). Let E G B(T). Then there exists E' £ A(N(T)) = A(N) 
such that \0\ = \E-E'\ < 1/T. 



Using Lemma 2.1, we obtain for all n with 1 < n < N, 



(39) \\T(n, 1;E+ ie)\\ < CN a exp(C(\f3\ + e)N a+1 ) = DN a , 

where D = Cexp(2C), since e = 1/T, \/3\ < 1/T and n < N = Similarly, 
for all n with — N < n < 0, 
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(40) \\T(n,0;E + ie)\\ < DN a . 



It follows from (B3h-(B4) and (B9h-(40) that 



(41) |0(n + 1)| 2 + |0H| 2 > KN- 2a (\d>(2)\ 2 + I0(1)| 2 ), 
for all E g B(T), 1 < n < N, and 

(42) \<j>{n + 1)| 2 + |0(«)| 2 > KN- 2a mi)\ 2 + |0(O)| 2 ), 

for all E € B(T), —N < n < 0. Here, > is some uniform constant. Now with 



l|)-(|42j), we can estimate, for any E £ B(T), 

N -N/2-1 

£ (|0(n + l)| 2 + |0(n)| 2 ) > E + E 

|n|>JV/2+l n=N/2+l n=-N 

(43) > B(A72-2)iV- 2Q (|0(O)| 2 + 2|0(l)| 2 + |0(2)| 2 ) 

> dn i - 2 <*(\<p(o)\ 2 +\m\ 2 + \m\ 2 ) 

with uniform constant D > 0. Here we have assumed that T is large enough, so 
that N > 8, for example. We now use equation ( |28| ) for the resolvent for n = 1: 

0(2) + 0(0) + — E — ie)0(l) = 1. 

Since £ g B(T),A(N) C [-if, if] and e = 1/T < 1, we get 

|0(2)| + |0(1)| + |0(O)|> ' 



|V(l)|+if + 2 
Together with ( [43] ) this gives 

(44) £ (|0(n+l)| 2 + |0(n)| 2 )> 7 iV 1 - 2Q , 

|n|>Af/2+l 

where 7 > is some uniform constant depending on |V(1)|, K, D. It is clear that 
the left-hand side of ( f44| ) is bounded from above by 

2 E i^)i 2 - 

|n|>W/2 

Thus, we obtain the following lower bound: For every E £ B(T), 
J2 \R(E + is)S l (n)\ 2 >lN 1 - 2a . 

\n\>N/2 

Therefore, 

(45) / V \R(E + ie)S 1 (n)\ 2 dE > f J TV 1 - 2 " dE = ^-\B{T)\N x - 2a . 



' |n|>JV/2 



Summation over {n : \n\ > N/2} in d3§| ) together with ( f45| ) proves (|l7|). Since 
|X|f(t) > (N/2Y \eM~m5i{n)\ 2 : 



\n\>N/2 
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the bound (18) follows immediately. 

In the half-line case, the proof is simpler. The equation (|3(]) for the resolvent 
for n = 1 gives the uniform lower bound for ||<j>(l)||, and using (|4l|), we prove the 
two statements of the Theorem. □ 



Proof of Corollary 1.1. One takes A(N) = {E } for all N. Since \B(T)\ = 2/T, 



the result follows directly from (|18|). □ 



Proof of Corollary U. One takes A(N) = [E - N'^^Eq + N' 1 ' 6 ]. One sees 



easily that fll6|) holds with a = 0. Therefore, N(T) = T and \B(T)\ > \A(N(T))\ = 
2T -i/e Thc bound Q yields the resu i t , □ 



Proof of Corollary \1.3[ . It is easy to see that the operator H satisfies the hypothesis 
of Theorem ]TJwith thc same K , a, A(N), and some appropriately adjusted constant 
C. This yields (pi) and (ph. □ 



3. The Fibonacci Model 

In this section we apply Theorem [l] to the Fibonacci model. We refer the reader 
to [Q for background information on this model and related ones. While Theorem]!] 
already gives a non-trivial dynamical bound when the set A(N) consists of a single 
point, and such an input is known for the Fibonacci model j22fl , we shall nevertheless 
identify the natural set A{N) of energies for which one can prove the local power- 
law bounds required by Theorem 2. On the one hand, this improves the dynamical 
bound, since we obtain a larger lower bound for the measure of B(T), and on the 
other hand it illustrates nicely how the sets, on which one can prove local power-law 
bounds (they will turn out to be spectra of periodic approximants) shrink down 
to a zero-measure set (the spectrum of the Fibonacci Hamiltonian) of energies for 
which the power-law bounds hold globally. 

Let us first present some basic notions and results for the Fibonacci potential 
given by (|2l|). The transfer matrices T(n, m; z) are defined as in (p6|), and we shall 
sometimes write T\(n, m; z) to make their dependence on the parameter A explicit. 

Define the sequence (Ffe)fc>o of Fibonacci numbers by 

F = 1, F 1 = 1, F k+l =F k + F fe _! for k > 1. 

Let 

x k {E,X) =tTM k (E,\), 

where 

M k (E, X)=T x (F k ,l;E). 

We will often leave the dependence of x k or M k on E, A implicit. The matrices M k 
obey the recursion 

(46) M fe = M k - 2 M k -i 
which yields 

(47) x k+ i = x k x k -i - x k - 2 
for their traces. This in turn gives the invariant 



(48) 



x l+i + x l+ x l-i - Xk+\x k x k -i = 4 + A 2 for every k e N. 
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The identities (^-(Q were proved by Siito in j35|. 
For fixed A, define 

a k = {EeR: \x k (E,X)\ < 2}. 
The set crj. is actually equal to the spectrum of the Schrodinger operator H whose 
potential 14 results from V in ( pi] ) by replacing a by Fk-i/F). (see pl|). Hence, 
14 is periodic, Ofc C M, and it consists of Fk bands (closed intervals). 

Lemma 3.1. For every E € o~k, we have 

(49) \xi(E, A) | < C x forO<i<k. 

Proof. It was shown in (35) that for every m G N, 

(T m U <7 m +l C (7 m _l U <7 m . 

Thus, if E £ crfc, then |xfc(-E, A)| < 2 and for every i < k, we have that either 
A)| < 2 or m&x{\xi~i(E, \)\,\xi+i(E, A)\} < 2. In the former case, the 
claimed bound clearly holds. In the latter case, use the invariant ( f48| ) to establish 
the claimed upper bound for \xi(E, A)|. □ 



Lemma 3^ is the key tool in identifying a natural set of energies for which the 
desired power-law bounds on transfer matrices up to a certain distance from the 
origin hold. 

Proposition 3.2. For every X, there is a constant C such that for every HN, 
every E £ a^, and every m with — Fk + 1 < m < (m 7^ 0), we have 

(50) \\T x (m,l;E)\\<C\m\ a , 
where 

log(C A (2C A + l) 2 ) 
51 a= — . 

log LO 1 

Proof. We first note that due to the symmetry V(—n) = V(n — 1), n > 2 of the 
potential |35), we can restrict our attention to 1 < m < Fk- 

Our principal strategy is to modify the approach of [^2| as to treat a larger set 
of energies, while proving bounds only in finite regions. Essentially, their proof of a 
power-law upper bound for energies in the spectrum can be turned into the claimed 
result by noting that for each m, boundedness of traces is only needed for j's with 
Fj < m, and the previous lemma established this fact for energies in o^. , where 
k = max-jj : Fj < m}. 

For the reader's convenience, let us be more concrete. Fix A, k, and E E at- By 
symmetry (V(— n) = V(n— I), n > 2; sec |3j|), we can restrict our attention to the 



case of positive m. The proof of (50) will be split into several steps. 
Step 1. For < i < fc, we have 

(52) \\Mi\\ < (Cx)\ 

This follows from 

Mi = M^ 2 M^ X = Mi^Xi-tl - Mr_\) = s i _ 1 M i _2 - M~i 3 
along with || Mr_\ || = ||Af 4 _ 3 ||. 
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Step 2. For i,j G N with i + j < k and i > 2, we have 

(53) MiM i+ j = P, (1) A/ l+ , + Pj 2) M i+J -_ x + pj 3) M i+j - 2 + pf ) I, 

where, for i = 1, 2, 3, 4, P. = P- a;,, . . . , Xi+j) is a polynomial of degree at 

most j, and we have 



(54) £|Pf < (2C A + 1)^, 
i=i 

where the polynomial |Pj | results from P- by replacing all the coefficients by 
their respective absolute values. 

To see this, given (ff9|), one can literally redo the proof of Lemma 5 in |2^] since 
it only uses the trace bound for indices bounded by k. 

Step 3. For 1 < m < Pfc, we have 

(55) \\T x (m,l;E)\\ <d m ", 
where 

d = C A (2C A + l) 2 

and rriN appears in the unique coding of m in terms of the Fibonacci numbers, 

N 

m = ^2 F "u > F m < F mi < ■ ■ < F mN , mi - mi-i > 2. 
;=o 

Clearly, = max{i : Pj < m} and hence mjv < k. 

The estimate (^5|) can be proved in the exact same way as in ]2^] , given Steps 1 
and 2 above. 

Step 4- Conclusion of the proof: The Fibonacci numbers P$ behave asymp- 
totically like uj~' 1 = [(1 + \/5)/2] 1 and in particular, there is a constant D\ such 
that 

Fi > D\U)~ l for every i > 1. 
This means that the number mjv defined in Step 3 obeys 

log 772 

TOat < ; r + £>2 

logw 1 

for some suitable constant D2, which in turn implies 



log' 



-D-2 



||T A (m,l;P)|| <d"^ 
and hence (pfj) for a constant C depending only on A and a as in (|5l|). □ 



We therefore have found a natural candidate that serves the purpose of the set 
A(N) in Theorem |l|. Namely, given some N, choose k such that Pfe_i < N < Fk 
and set A(N) = a k . 

Next, we prove a lower bound for the Lebesgue measure of ak- More precisely, 
we obtain lower bounds for each one of the Fk intervals it is made of. The results 
presented here are essentially contained in Raymond |3l| . We give a somewhat 
more streamlined presentation in the spirit of pfl . 

From now on, we shall always assume 
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(56) 



A > 4, 



since we will make critical use of the fact that in this case, it follows from the 
invariant ( [48| ) that three consecutive traces cannot be bounded in absolute value 
by 2: 



(57) 



VA >4,VF,fc : max.{\x k (E,\)l\x k+1 (E,\)l\x k+2 (E,\)\} > 2. 



Following 1 26 , we call a band If. C Cfc a type A band if I k C cr k —i (and hence 
Ik H (crfe + i U er/c-2) = 0). We call a band If. Cat a type B band if Ik C cr/t-2 (and 
therefore Ik H CTfc_i = 0). 

From (|57j), one gets the following (Lemma 5.3 of j26fl , essentially Lemma 6.1 of 

0): 



Lemma 3.3. For every A > 4 and every k G N. 

(a) Every type A band Ik C <7fc contains exactly one type B band Ik+2 C 0^+2, 
and no other bands from o~k+i, o~fc+2- 

(b) Every type B band Ik C o~k contains exactly one type A band Ik+i C (Jk+i and 
two type B bands from Cfc+2, positioned around Ik+i- 

We will also need the following lemma (Lemma 5.4 of essentially Proposi- 
tion A. 2 of |3l)). 

Lemma 3.4. Let the functions f±(x,y,X) be defined by 

1 



f± (x,y,X) = - [xy ± V4A 2 + (4-o: 2 )(4-y2) 



For A > 4 and —2 < x, y < 2, we Ziawe 



<9x 



(a,y,A) 



dy 



(x,y,X) 



< 1. 



Equipped with the previous two lemmas, we are now in position to prove the 
following: 

Lemma 3.5. For every A > 4, the following holds: 

(a) Given any (type A) band Ik+i C o~k+i lying in the band Ik C o~u> we have for 
every E G h+i, 

x' k+ x(E) 



x'k(E) 



< A + 11. 



(b) Given any (type B) band Ik+2 C o~k+2 lying in the band Ik C o~k> we have for 
every E G h+2, 



< 2(A + 11) 



<(E) 

Proof, (a) Differentiating (|47]) and dividing by x' k , we get 



(58) 



b k+l 



Xk-1 



•KkXfc_ 



L k-2 



Using the invariant (48) and Lemma 



3.4. 



, we obtain for every E G Jfc+i, 
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\Xk- 



\f±{Xk,Xk-2,X)\ 
1 



x k x k - 2 ± a/4A 2 + (4 - 4) (4 - x\_ 



< 



i ^4+ V4A2 



64 



< A + 6. 



Lemma 3.4 implies 



— (x k ,X k ^ 2 , X)x' k + -^-(x k ,X k -2, X)x' k _ 2 



(59) 

From Lemma 5.5 of [M, we infer, again for E G Ik+i> 



< 141 + 14-2 1 



(60) 

Using ©, ©, and ©, we get 



°fc-2 



< 1. 



< A + 6+|x fe |- 



"fc-11 



"k-2\ 



"k-2\ 



< A + 6 + 2^1 

< A + 11. 

(b) We consider two cases. Assume first I k+ 2 n er^-i = (and so I k+ 2 C <r k -2)- 
Given 

X 'k+2 = x k+l x ' k + x' k+1 X k — X k _ x 

and 

/ _ 1 1 / / 

^fe+l — x kX k _ 1 + X k X k -i — X k _ 2 i 

we find 



J k+2 



2x k+ i ~ X k ^2 + (x k - 1) — 7 



k-1 



x k - 



J k-2 



Similarly to the previous argument, 

\xk+i\ < A + 6 



and 

which yields 



"fc+2 



l4-il < I4I + I4- 2 I, 

14-2 I 



"k-2\ 



< 2A+ 12 + 2 + 3^1 

< 2A + 22. 

Let us now assume I k +2 C o k -\ (and so H o k -2 = 0). We proceed analogously 
and obtain 



X k+2 . X k+1 
— = X k+ i + X k — 



v k-l 
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and 

which then yields 



\%k+i \ < A + 6, 

l4+i I < 141 + l4-i I. 



b k+2 



Xu + x 



fc-il 



^fe-il 



< A + 6 + 2- 

\ x 'k\ \ x k\ 

< A + 11, 

concluding the proof. □ 
This yields, as a very rough estimate, that for A > 4 and E G ak, we have 

(61) I4(£,A)| < C(2A + 22) fc , 

which can easily be turned into a lower bound on bandwidths, as the following 
proposition shows. 

Proposition 3.6. For every A > 4, the set ak consists of F k disjoint closed inter- 
vals, each of which has Lebesgue measure bounded from below by 



(62) 

In particular, we obtain 
(63) 

where 7 is as in (pJ). 



> 



C(2A + 22)* 



4 

C' 



Proof. Since o~k is the spectrum of a periodic Schrodinger operator with period 
Fk, the first statement is immediate. That the gaps are open was shown in fjlf . 
The estimate ( |62| ) for the measure of one of these Fk bands can be seen as follows 
(cf. Jul): On such a band Ik, Xk(E, A) runs monotonically from ±2 to =f2. Hence, 
by ©, 

' \x' k (t,X)\dt< \I k \-C(2X + 22) k . 



The estimate (|63| ) then follows from this and the exponential behavior of the se- 
quence (F fc ) feeN . □ 

We have established the input to Theorem [I]. Note that in any event, we have 
\B(T)\ > 2/T. In the case where A > 4, Proposition |3.6| improves this factor in 
( fl7|) and (|l8|). We therefore proceed with the 

Proof of Theorem [|. (a) This follows from Proposition |3.2| and Corollary |1.1| . 

(b) Given T, we let as above N(T) = T 1 ^ 1 ^ and we choose k such that F k ~i < 

N(T) < F k . We let A(N(T)) = a k . Thus, from Proposition O we get 



\B(T)\ > ±FP > 1(2^.! 



4 4 — T 

)-7 > —2^N(Ty = — 2~ 1 T~ 

c c 



Now d|) yields 



l*D(r)>C2-^r^ 
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for T large enough and hence (3 (p) > (p — 7 — 3a)/ (1 + a). □ 
4. The Period Doubling Model 

In this section we investigate the period doubling model and, in particular, prove 
Theorem [| This will be done by first proving linear bounds on transfer matrix 



1.1 



norms for a finite/countable set of energies and then alluding to Corollary 

On the alphabet A = {0, 1}, consider the period doubling substitution 5(0) = 01, 
S(l) = 00. Iterating on 0, we obtain a one-sided sequence u = 01000101 . . . 
which is invariant under the substitution process. Define the associated subshift 
ilpd to be the set of all two-sided sequences over A which have all their finite 
subwords occurring in u. For A £ R and u> <E £l p d, we define the potential V\ tU by 
Vx, u (n) = Xuj„. 

If w = w% . . . Wi with Wi G {0, 1}, we define 

T x {w;E) = A x (wi;E) x ••• x A x { Wl ;E), 

where for a G {0, 1}, 

a 1 r-.\ f E — Xa —1 
A x (a;E)=l 1 Q 

We let 

T jfi) =T <f)(E,\)=T x (S k (0);E), T, (1) - T^(E, A) = T x (S k (l):E). 
In the following, we will leave the dependence of K , on E and A implicit. 



We also define 



x k = trT fc (0) , y k = trT^. 



It follows from the substitution rule (and is easy to check) that 



T (o) _ T (iW0) T (i) _ T (o) T (o) 

1 k+l — 1 k ± k 1 1 k+l — 1 k ± k 



and 



(64) = XkUk - 2, yk+i = x 2 k - 2. 

The relation (^i|) is called the period doubling trace map. 

By virtue of Corollary 1.1, Theorem || follows once we find an energy Eq such 
that ( |l9| ) holds with a — 1. We establish this first and then discuss later what the 
natural set of energies is for which one can establish a bound like (|l9|) . 

Proof of Theorem ^. Clearly, we have 



T (o) _ 



E -1 
1 



and 



T (i) 



E-X -1 
1 



Notice that for Eq =0, we have 



T (o) _ ( -1 \ (i) _( -X -I 
J ° - \ 1 I ' ~ I 1 



Thus, for this choice of the energy, we get 
T (o) _ T (i) T (0) _ 



X 














1 


-i 


J 


j 


-( 





-i 
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and 



T (i) _ T (o) T (o) _ 
1 i — 1 o 1 o — 


















1 






/ 


■( 





-J 



In particular, Tj and Tj commute. Notice that 



(if)' 



1 -nA 
1 



Every subword roofu can be partitioned into a product of blocks of the form 5(a) 
or 5(6), up to a possible prefix/sufhx of length one. We therefore get, for every 

LU G Slpd, 



\\T(n,m;E )\\<Cl[V2 + ~\n- 



where 



C 



A 







(1 





2 



< V2 + A. 



We can now apply Corollary 1.1. □ 



We now study the set of exceptional energies (i.e., where we get commuting 
transfer matrices at some level) a little further. We prove the following: 

Lemma 4.1. Fix some coupling constant A and some k G No- If E is such that 
x n = 0, then there is a constant a\,k such that 



-1 



(65) ^fc+i * s conjugate to 

and 



(66) r«=(-J _° 

Proof. If a;*.. = 0, then ( |64| ) gives xt+i — —2. Hence, ( |65| ) is immediate. Moreover, 
( |66| ) follows from the Cayley-Hamilton theorem since we have 

rp{l) _ rp{0) T {Q) _ T (0) T , _ f -1 \ 

J fc+1 ~~ k ± k — X k--L k I Q — X / ' 

□ 

This allows us to prove: 

Proposition 4.2. Fix some coupling constant A and some fc G No- For every roo< 
£■0 of Xk, we have 

\\T(n,m;E )\\<Cl k (V2+^\n-m\), 
where C\_k is some suitable constant. 

Proof. The argument is virtually the same as the one used in the proof of Theo- 
rem [| □ 

Notice that Xk is, as a function of E, a polynomial of degree 2 k which has exactly 
2 k roots. In particular, there is a countable set of energies where linear bounds on 
the growth of transfer matrix norms can be proved. 
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5. The Thue-Morse Model 

In this section we investigate the Thue-Morse model and, in particular, prove 
Theorem || This will be done by exhibiting a set of energies for which the transfer 
matrix norms are bounded and then alluding to Corollary |l.l| . 

On the alphabet A — {0, 1}, consider the Thue-Morse substitution 5(0) = 01, 
5(1) = 10. Iterating on 0, we obtain a one-sided sequence u — 01101001 . . . which 
is invariant under the substitution process. Define the associated subshift O tm to 
be the set of all two-sided sequences over A which have all their finite subwords 
occurring in u. For A £ R and u> € ^tmi we define as in the period doubling case 
the potential Vx >u) by V\. w (n) = Xuj n . 

We can now define transfer matrices in the same way as above, that is, 

Tfc (0) =T ^(E,X)=T x (S k (0);E), T« = T« (E, A) = T A (5 fc (l); E). 
Again, we will leave the dependence of on E and A implicit. We define 

x k = trT fc (0) , y k = trT fc (1) . 
It is clear that x k — yk for fc > 1 and it follows from the substitution rule that 



T (o) _ T (i) T (0) T (i) _ T (o) T (i) 

1 k+l — 1 k ± k ^ ± k+l^ 1 k ± k 



and 



(67) x k +i = x|_x(x fc - 2) + 2 for k > 2. 

The relation ( |S7| ) is called the Thue-Morse trace map. 

By virtue of Corollary LI , Theorem [| follows once we find an energy Eq such 



that (H) holds with a = 0. Let £ k = {E : x k = 2} 
Proposition 5.1. If k > 3 and E € £k\^2> then 

(68) 2f>=lf>=(j J 

Proof. This can be extracted from [Q. For the reader's convenience, we give a short 
proof of this fact. Notice that it follows from (|67]) that for k > 3, we have Xk = 2 
if and only if Xk-i = 2 or Xk-2 = 0. Iterating this, we get that Xk = 2 holds if and 
only if X2 — 2 or xj = for some 1 < j < k — 2. Thus, if fc > 3 and E £ £k\£2, 
then Xj = for some 1 < j < k — 2. Using this and the Cayley-Hamilton theorem, 
we obtain 



T (0) _ (0)^(1)^(1)^(0) 
1 j+2 ~ 2 j 1 j 1 3 2 j 



(0) 
3 



= _ r (o) T (o) 
i j 

- -{*flf- 

= I 

and, similarly, T^ 2 = I, This yields @. □ 
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Proof of Theorem^ Given Proposition 5.1, the claim follows as in the proof of 



Theorem || from Corollary □ 

For example, the reader may verify that if E is chosen such that E(E — A) = 2, 
then ( p8| ) holds for k — 3 by a straightforward calculation. This observation already 
suffices for an application of Corollary [P| . 

We conclude with a remark about the special energies exhibited by Proposi- 
tion 5.1. It follows from (p7\) that for k > 2, Ek C £k+i- Axel and Peyriere 



show that the the union of the sets £k is dense in the spectrum of the Thue-Morse 
Hamiltonian IpJ. 
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